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Abetrect 

An  embedded  mesh  system  for  transonic  flow  analysis  is  described.  The 
method  has  been  applied  to  the  finite-volume  potential  flow  calculation  of 
flows  around  wlng-body-tall  configurations.  A  local  C-tjrpe  mesh  is  fitted 
to  a  horizontal  tall  mounted  on  a  wing-body  around  which  a  wraparound  C-type 
mesh  has  also  been  used.  Sample  calculations  for  a  modified  A-7  configura¬ 
tion  and  another  wing-body-tall  configuration  are  given.  The  results  show 
rapid  convergence  and  excellent  details  of  pressure  distribution  over  the 
entire  configuration  in  both  cases. 
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1.  Introduction 

The  finite-volume  algorithm  for  calculating  transonic  potential  flow 
around  an  arbitrary  geometry  was  orglnally  developed  by  Jameson  and  Caughey 
(1977).  This  method  Is  now  used  as  the  base  for  computational  codes  that 
analyze  the  flow  around  realistic  aircraft  configurations.  One  of  the  most 
popular  of  these  codes  Is  the  FLO-30  code,  originally  developed  by 
Flow  Research  Company  In  cooperation  with  Professors  Jameson  and  Caughey 
(Mercer  et  al.,  1980),  which  uses  a  C-type  mesh  system.  The  finite-volume 
algorithm  is  so  flexible,  however,  that  other  equally  effective  and  more 
flexible  codes  can  be  constructed  using  other  types  of  mesh  systems.  For 
example,  the  WBT  code  has  been  developed  by  Flow  Research  Company  for 
Arnold  Engineering  Development  Center  for  simulating  transonic  flow  around 
a  wing-body  combination  in  a  wind  tunnel.  This  code  has  been  used  to  model 
the  effects  of  arbitrary  wall  boundary  conditions  on  wind  tunnel  models. 
Good  agreement  can  be  seen  between  results  from  the  WBT  code,  using  a 
constant  pressure  condition  at  the  distant  wall,  and  FLO-30.  This  agree¬ 
ment  verifies  that  the  finite-volume  algorithm  can  be  used  on  a  wide 
variety  of  grid  systems. 

It  is  now  possible  to  add  other  components  of  an  aircraft  to  these 
basic  wing-body  codes  to  enhance  their  capabilities.  In  the  present  work, 
a  horizontal  tall  is- incorporated  in  the  basic  configuration.  Horizontal 
tails  with  substantial  lifting  force  can  couple  fairly  strongly  with  the 
flow  around  the  wing,  particularly  when  the  wing  and  tail  are  closely 
coupled  as  in  many  military  aircraft. 

In  this  report,  the  mesh  generation  scheme  for  our  wing-body-tall 
combination  is  described.  We  show  that  if  the  C-type  mesh  is  required 
around  the  horizontal  tall  to  resolve  the  leading-edge  pressure  detail,  it 
is  topologically  impossible  to  generate  a  unified  mesh  system  that  includes 
both  wing  and  tall;  an  embedded  mesh  must  be  used  for  the  horizontal  tall. 
The  solution  is  then  obtained  through  interactive  iterations  between  the 
main  wing-body  calculations  and  the  tall  calculations.  Examples  of  these 
calculations  are  presented. 
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2.  Mesh  Generation 

The  existing  wing-body  codes,  such  as  FLO-30  and  WBT,  use  a  C-type  mesh 
system  around  the  wing.  The  wraparound  mesh  lines  are  carried  to  the  down¬ 
stream  infinity.  If  a  wraparound  C-type  mesh  is  required  for  the  tall 
region,  these  two  families  of  C-type  curves  will  be  compatible  downstream 
of  the  tall  but  Incompatible  in  the  space  behind  the  wing  and  in  front  of 
the  tall.  This  means  that  a  single  mesh  that  simultaneously  wraps  around 
the  wing  and  the  tall  cannot  be  generated.  Therefore,  two  mesh  systems  must 
be  employed,  one  for  the  wing  and  one  for  the  horizontal  tall,  and  the  tall 
system  must  be  embedded  into  the  wing  system.  As  a  means  of  defining  the 
outer  boundary  of  the  tall  region,  a  set  of  grid  lines  from  the  wing  grid 
system  is  used  as  shown  in  Figure  1.  This  procedure  creates  two  computa¬ 
tional  modules:  a  wing-body  calculation  and  a  tall  calculation.  The 
boundary  of  the  wing-body  calculation  extends  into  the  interior  of  the  tail 
region.  This  overlapping  of  the  two  grids  allows  for  the  transfer  of 
information  between  the  two  calculations  as  discussed  in  Section  3. 

2.1  Wing-Body  Mesh 

The  WBT  code  developed  for  a  wing-body  combination  in  a  wind  tunnel  was 
chosen  in  favor  of  the  FLO-30  code  as  the  basic  code  for  incorporation  of  a 
horizontal  tall.  There  are  several  reasons  for  this  choice.  Tbe  FLO-30 
code  uses  a  cylindrical-type  mesh  which  has  ’’spanwise"  coordinate  surfaces 
wrapped  around  the  fuselage.  This  system  makes  it  very  difficult  to  model 
tall  configurations  which  lie  near  or  off  the  edge  of  the  fuselage,  such  as 
is  often  the  case.  Also,  the  vortex  wake  trajectory  from  the  wing  becomes 
difficult  to  model  with  the  cylindrical-type  mesh.  This  is  due  to  the  fact 
that  the  wake  must  lie  on  a  grid  surface  and  so  must  be  displaced  to  a 
boundary  of  the  tall  mesh.  In  general,  the  wake  cannot  pass  through  the 
tall  since  there  is  no  mesh  line  in  a  C-type  mesh  that  goes  from  the  upstream 
boundary  to  the  downstream  boundary  (except  the  one  which  passes  through  the 
tall) .  If  the  vortex  wake  is  distorted  to  the  tail  region  boundary  on  a 
cylindrical-type  mesh,  the  wake  will  not  only  be  displaced  vertically  but 
also  horizontally.  The  horizontal  displacement  can  produce  large  errors  for 
the  downwash  effect  on  the  horizontal  tall.  The  WBT  code  uses  a  mesh  system 
that  results  from  transforming  the  fuselage  cross  sections  to  a  vertical 
slit.  The  horizontal  tall  is  easily  modeled  as  a  tall  on  a  wall  in  the 
transformed  space  whether  it  is  on  the  fuselage  or  not.  Also  a  spanwise 
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station  is  roughly  a  spanvise  slice.  The  displacement  of  the  wing  vortex 
wake  in  the  spanwise  plane  moves  the  vortex  wake  upward  without  signifi¬ 
cantly  changing  the  spanwise  location.  Lastly,  a  code  that  can  be  used  for 
t>oth  free— air  calculations  and  wind  tunnel  simulations  will  be  extremely 
valuable  both  as  a  design  tool  and  as  a  tool  to  assess  wind  tunnel  test 
results. 

The  procedure  for  generating  a  wing-body  mesh  in  the  WBT  code  can  be 
summarized  as  follows: 

(1)  The  first  step  is  to  transform  a  series  of  circular  cross 
sections,  with  centers  and  radii  determined  by  the  height  and 
location  of  the  fuselage  cross  sections,  to  vertical  slits  by  use 
of  the  well-known  Joukowsky  transformation.  Under  this 
transformation,  a  noncircular  fuselage  will  become  a  distorted 
slit. 

(2)  Next,  a  shearing  transformation  in  the  spanwise  direction  is  used 
to  remove  this  distortion  of  the  fuselage. 

(3)  The  computational  domain  in  this  Joukowsky-shear  transformed  space 
is  then  sliced  into  spanwise  sheets  on  which  the  wraparound  C-type 
mesh  will  be  generated. 

(4)  On  the  spanwise  station  wrapping  around  the  fuselage,  two 
boundaries  must  be  placed  on  mesh  lines:  one  is  the  intersection 
of  the  fuselage  with  the  wing  root,  and  the  other  is  the  upper 
and  lower  edges  of  the  fuselage  that  has  been  transformed  to  a 
slit.  Both  of  these  boundaries  must  be  mesh  lines  of  the 
wraparound  family.  To  be  consistent  with  the  topology  of  the 
fuselage  station  just  described,  images  of  the  "fuselage"  are 
defined  for  the  other  spanwise  stations.  These  images  are  defined 
so  that  the  mesh  system  continuously  varies  from  the  fuselage 
shape  at  the  Innermost  station  to  an  open-elliptical  shape  at  the 
tip  station  and  beyond. 

(5)  For  wind  tunnel  applications,  an  additional  boundary  is  required. 
This  is  a  rectangular  boundary  at  the  outer  limit  of  the  computa¬ 
tional  domain  corresponding  to  the  tunnel  wall.  For  free-air 
calculations,  this  boundary  is  assumed  to  be  a  large  distance  away 
from  the  configuration. 


( 
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It  is  not  possible  to  generate  a  family  of  C-type  curves. 

Including  three  known  curves,  with  a  conformal  transformation. 
Thus,  the  application  of  the  mesh  generation  scheme  using  an 
elliptic  equation  is  ruled  out  when  it  is  cot.pared  to  the  follow¬ 
ing  simple,  algebraic,  mesh  generation  scheme. 

A  family  of  "superellipses'*  is  generated  by  the  following  formula: 


m 

m 

(-)  + 

\a  / 

\b  / 

where  the  parameter  a  is  the  semi-major  axis,  b  is  the  semi- 
minor  axis,  and  m  controls  the  curvature  of  the  superellipse. 

For  example,  as  m  approaches  infinity,  the  superellipse 
approaches  a  rectangle.  The  family  of  superellipses  used  in  this 
work  is  characterized  by  the  following  properties: 

(a)  The  parameter  b  varies  from  the  half-thickness  of  the 
wing  at  the  innermost  curve,  to  the  maximum  half-height 
of  the  "fuselage"  image  at  a  designated  mesh  index,  and 
to  the  "tunnel"  half-height  at  the  outer  boundary. 

(b)  The  semi-major  axis,  a  ,  varies  from  the  chord  length  of 
the  wing,  to  the  distance  between  the  nose  of  the 
fuselage  and  the  trailing  edge  of  the  wing  at  a 
designated  mesh  index,  and  to  the  distance  between  the 
upstream  boundary  and  the  trailing-edge  boundary. 

(c)  The  parameter  m  varies  from  2  at  the  innermost  ellipse 
to  infinity  at  the  outer  boundary. 

(d)  The  trailing  edge  of  the  wing  is  defined  as  the  origin 
of  the  family  of  superellipses. 

This  process  generates  a  family  of  curves  of  which  three  members 
have  approximate  thickness  ratios  of  the  wing,  the  fuselage,  and 
the  tunnel  wall  (outer  boundary),  respectively.  These  C-type 
curves  are  continued  horizontally  from  the  trailing  edge  to  the 
downstream  boundary. 


(7)  To  form  a  computational  mesh,  these  C-type  curves  must  be 

intersected  by  another  family  of  curves.  The  innermost  curve  is 
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divided  into  segments  so  that  the  segments  are  small  near  the 
leading  edge.  The  outer  boundary  is  also  divided  into  an  equal 
number  of  segments.  Starting  from  the  upstream  point,  the 
corresponding  points  on  these  two  curves  are  connected  by  straight 
lines.  The  intersections  between  the  straight  lines  and  the 
C-type  curves  are  computed.  These  intersections  define  a  mesh 
network. 

(8)  The  next  task  in  the  mesh  generation  is  to  make  it  a  body-fitted 
mesh.  This  is  accomplished  by  shearing  the  mesh  vertically  so 
that  the  Innermost  C-type  curve  is  moved  to  the  actual  wing 
profile  and  the  C-type  curve  designated  as  the  fuselage  mesh  line 
is  moved  to  the  actual  fuselage  image. 

(9)  Finally,  these  two-dimensional  meshes  are  transformed  to  the 
physical  space  by  an  inverse  Joukowsky-shearing  transformation. 

This  procedure  for  generating  a  wing-body  mesh  is  demonstrated  schematically 
in  Figure  2. 

This  mesh  generation  scheme  has  a  number  of  advantages  over  the  FLO-30 
scheme.  First  of  all,  the  scheme  can  be  used  for  free-air  calculations  as 
well  as  wind  tunnel  simulations,  when  the  latter  are  required.  Second,  the 
extension  of  the  mesh  upstream  is  an  input  parameter,  which  Insures  that  the 
upstream  boundary  conditions  are  applied  at  a  distance  sufficiently 
upstream.  It  has  been  observed  that  the  FLO-30  mesh,  particularly  that 
portion  near  the  futelage,  does  not  extend  very  far  upstream.  This  has 
caused  some  problems  in  the  accuracy  of  the  application  of  the  far-fleld 
free-stream  conditions. 

2.2  Tail  Mesh 

In  the  Joukowsky-sheared  space  of  the  wing-body  mesh,  a  region  in  the 
neighborhood  of  the  tail  is  cut  out.  This  region  is  defined  by  a  vertical 
line  in  the  wing-body  mesh  a  certain  distance  upstream  of  the  leading  edge 
of  the  tall  and  by  two  horizontal  mesh  lines.  The  upper  and  lower 
boundaries  of  the  tail  region  must  be  defined  sufficiently  far  away  from  the 
tail  so  that  the  supersonic  bubble  generated  by  the  tail  does  not  cross  the 
boundaries.  Due  to  the  disparity  in  the  mesh  spaclngs  between  the  wing-body 
mesh  and  the  tall  mesh  on  these  boundaries.  Interpolation  of  the  velocity 
potential  there  can  cause  errors  in  the  mass  conservation  formula.  This  can 


•  Use  Super-Ellipses 


Figure  2.  Mesh  Generation  Procedure 
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be  a  severe  problem  In  the  supersonic  region.  The  present  choice  of  the 
outer  boundary  above  and  below  the  fuselage  Image  will  create  some  "fin" 
effects  (see  Mercer  et  al.,  1980).  However,  the  tall  mesh  is  sufficiently 
dense  so  that  the  fin  effects  are  not  pronounced. 

The  tall  mesh  is  generated  In  essentially  the  same  way  as  the  wing-body 
mesh,  except  that  the  fuselage  Image  Is  ignored  and  the  outer  boundary  is 
that  given  by  the  wing-body  mesh. 

For  a  modified  A-7  airplane,  the  side  view  of  the  combined  wing-body 
and  tail  meshes  on  the  "fuselage"  spanwlse  station  is  shown  in  Figure  1. 

For  the  purpose  of  clarity,  the  wing-body  mesh  outside  of  the  fuselage  image 
is  not  plotted.  The  fuselage  in  the  tail  region  shows  some  small  "fins"  on 
its  top  and  bottom. 
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3.  Numerical  Calculations 

The  calculation  of  the  velocity  potential  Is  carried  out  by  using  the 
finite-volume  method  with  the  following  strategy  of  Interactive  computations: 

(1)  The  wing-body  calculation  Is  performed  with  a  region  of  the  tail 
cut  out.  The  boundary  points  at  the  Interface  extend  one  mesh 
line  Into  the  tall  domain.  The  potential  values  on  the  boundary 
points  are  obtained  by  Interpolation  from  the  Interior  points  of 
the  tall  computations. 

(2)  During  the  wing-body  calculations,  the  circulation  around  the  wing 
section  continuously  Increases  with  each  iteration.  The  change  in 
circulation  after  each  iteration  is  divided  between  the  upper  and 
lower  boundaries  around  the  tall  region  and  is  added  or  subtracted 
from  the  potential  values  on  the  boundary  points. 

(3)  After  a  few  iterations,  the  computational  mode  is  switched  to  the 

tail  calculation.  The  potential  at  the  boundary  points  in  the  interior 
of  the  wing-body  domain  Is  obtained  by  interpolation  from  the  wing-body 
calculations,  with  the  proper  potential  jump  subtracted  from  the 
potential  on  those  points  located  above  the  vortex  wake.  This  opera¬ 
tion  equivalently  places  the  wing  vortex  wake  on  the  upper  boundary  of 
the  tall  region.  The  downwash  and  the  averaged  spanwlse  flow  Induced 
by  the  wing  are  fairly  well  approximated  by  this  method.  Only  the 
velocity  shear  on  the  vortex  wake  Is  at  the  displaced  location. 

(4)  After  a  few  Iterations  in  the  tail  computation,  the  potential  on 
the  boundary  points  of  the  wing-body  mesh  Is  obtained  through 
Interpolation  with  the  proper  potential  jump  added. 

(5)  The  process  is  repeated  until  convergence. 

The  code  is  written  for  use  in  a  computer  with  virtual  memory  or  a  large-core 
computer  such  as  the  Cray  1. 

A  general  "field-point  Interrogator"  is  programmed  to  Interpolate  the 
velocity  potential  at  the  boundary.  This  interrogator  essentially  treats 
the  transformation  from  the  physical  space  (x,  y,  z)  to  the  transformed 
space  (1,  j,  k)  as  a  nonlinear  function.  Given  a  set  of  arbitrary  physical 
coordinates  (x^,  yQ,  Zq),  one  can  treat  the  (1,  j,  k)  space  as  a  continuous 
space  and  locate  a  noninteger  location  (1^,  j^,  k^)  in  the  computational 
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space  by  use  of  Newton's  method.  The  velocity  potential  at  that  point  can 
be  obtained  from  that  at  the  surrounding  mesh  points  using  the  trlllnear 
formula.  This  Interrogator  makes  the  transfer  of  Information  at  the  inter- 
facial  boundary  extremely  simple. 
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4.  Results 

The  code  described  In  the  previous  sections  has  been  applied  to  a 
tnodlfled  A-7  aircraft  with  a  high  wing  position.  The  tail  Is  modeled  as  a 
7.S-percent  symmetric  NACA  00075  airfoil  with  a  -6  degree  angle  of  Incidence 
relative  to  the  defining  coordinate  axis.  The  planforms  of  the  wing  and  the 
tail  with  their  frontal  views  are  given  In  Figure  3. 

The  computed  pressure  distributions  on  the  wing  and  the  tail  are 
displayed  in  Figures  4  through  12.  For  a  free-stream  Mach  number  of  0.85 
and  an  angle  of  attack  at  4.68  degrees,  the  pressure  distribution  on  the 
tail  does  not  show  any  shock  wave  effects.  Near  the  trailing  edge  of  some 
of  the  wing  sections,  however,  a  dip  in  the  pressure  plots  can  be  seen. 
Careful  examination  of  the  given  data  shows  that  a  concave  kink  Is  present 
on  the  upper  surface  of  the  airfoil  near  the  trailing  edge.  No  attempt  has 
been  made  to  smooth  the  kink  and  possibly  improve  the  appearance  of  the 
pressure  plot. 

The  code  was  also  applied  to  a  wing-body  geometry  given  by  SAAB  of 
Sweden.  The  configuration  has  a  high  wing  geometry  with  a  cylindrical  body, 
and  a  horizontal  tall  has  been  added.  The  tall  Is  a  NACA  0005  airfoil  at  a 
zero  degree  angle  relative  to  the  defining  coordinate.  The  planforms  of  the 
wing  and  the  tail  together  with  their  frontal  views  are  given  in  Figure  13. 
The  pressure  distributions  on  the  wing  and  the  tall  are  given  in  Figures  14 
through  22. 

The  iterative  process  In  the  code  Is  organized  into  outer  wing-tail 
Interactive  iteration  loops  and  Inner  iterations.  We  found  that  good 
convergence  can  be  achieved  with  five  Inner  Iterations  for  both  wing-body 
calculations  and  tail  calculations.  Since  the  convergence  Is  fairly  rapid, 
the  total  number  of  Iterations,  l.c.,  the  number  of  outer  Iterations  multi¬ 
plied  by  the  number  of  inner  Iterations,  is  comparable  to  the  wing-body 
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g  and  Tail  Pianforma  for  A-7 


Flow  Research  Report  No.  227 
April  1982 


s  » 


*»(«***»< 


xxxxxxM'^ 


«  •  K  + 


M 

-I-  K 


X  X 


4 


4 

4- 


4  * 


K 

* 


fi-7  TRSON  HING  N2  H/B  FUS  FIB  WING  AT  T=120. 
MflCH  0.  850  RLPHR  4.  680 

Z  29.02  CL  0.1793  CD  0.0418 


Figur*  4.  A-7  Wing  Prastur* 
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R-7  TRSON  HING  H2  H/B  FUS  FIB  WING  AT  T=120. 
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Figures.  A-7  Wing  Pressure 
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Figur*  8.  A-7  Wing  PrMsur* 
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fl-7  TRSON  WING  W2  N/B  PUS  FIB  WING  RT  Y=120. 
MRCH  0. 850  RLFHfl  4.  680 
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Hgur*  10.  A-7  Tail  Praatura 


Flow  Research  Report  No.  227 
April  1982 


-20- 


+ 

+ 


-*■ 

+ 

4. 


fl-7  TR30N  HING  W2  W/B  FUS  FIB  WING  RT  ^=120. 
MflCH  0.650  RLPHfi  4.680 

Z  74.16  CL  -0.2295  CD  -0.0105 


Figur*  1 1 .  A-7  Tail  Prassure 
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5.  Conclusions 

The  present  work  demonstrates  that  the  basic  wing-body  code  can  be 
extended  to  Include  other  major  components  of  a  realistic  aircraft  by  using 
an  embedded  mesh  with  interactive  computations.  The  interaction  between 
modules  does  not  substantially  affect  the  rate  of  convergence.  Results  for 
closely  coupled  wing-body-tail  configurations  show  that  the  method  is  stable 
and  produces  the  expected  trends.  Comparison  with  wind  tunnel  data  for  wing- 
tail  configurations  should  be  made  when  available  to  test  the  accuracy  of  the 
method.  The  program  can  provide  designers  with  an  accurate  method  for 
assessing  tall  designs  under  transonic  flight  speeds.  For  very  closely 
coupled  configurations,  however,  the  method  may  fall  because  there  may  not 
be  enough  space  between  the  wing  and  the  tail  for  generating  an  embedded 
C-type  mesh  around  the  tall.  Another  mesh  system,  such  as  a  unified  grid 
with  an  H-type  mesh,  may  be  required  (see  Huynh  and  Jou,  1982).  However, 
the  concepts  of  overlapping  meshes  and  interactive  iterations  demonstrated 
here  can  be  of  great  value  in  any  grid  system  where  local  accuracy  is 
required. 
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